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- ^The  task  of  finding  shear  distributions  in  damaged  composites  analytically 

by  means  of  the  theory  of  elasticity  is  prohibitively  difficult,  so 
approximate  techniques  are  generally  employed.  The  most  widely  used  is 
shear  lag  theory.  This  theory  contains  a  parameter  representing  the  shear 
coupling  between  fibers that  has  been  difficult  to  evaluate  in  the  past. 

The  present  paper  evaluates  this  parameter  for  a  linear  array  of  fibers. 

In  the  limiting  case  of  small  fiber  diameter  to  fiber  spacing  ratio 
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20.  Abstract 


and  large  matrix  thickness  to  fiber  diameter  ratio  an  analytical  solution 
is  possible  and  is  provided  in  this  paper.  For  other  geometries  the 
solution  is  obtained  experimentally.  It  is  found  that  the  assumption 
usually  made  for  the  parameter  is  rather  far  off  for  most  geometries  of 
practical  interest.- — 
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distance  between  fiber  or  rod  centers 

distance  between  nearest  fiber  pair  in  a  rectangular  array 

h/ d  for  the  interaction  between  nearest  fiber  pair  in  a  rectangu¬ 
lar  array 

radius  of  a  fiber  or  a  rod 

lamina  thickness 

axial  displacement  of  the  i 9 th  fiber 

axial  coordinate  t  Hi-.- 


load  transfer  to  the  i 9 th  fiber  per  unit  length 
axial  line  force  per  unit  length 


tm 


shear  modulus  of  the  matrix 

effective  shear  stiffness  of  the  composite 

current  generated  by  the  power  supply 

a  large  distance  approaching  infinity  at  which 
the  field  point  is  essentially  zero 


resistance 


a 


ratio  of  the  distance  between  next- to-nearest  fiber  pair  to  the 


distance  between  nearest  fiber  pair  in  a  rectangula 


o  resistivity  of  the  electrolyte 
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p  electrical  potential 


INTRODUCTION 


Studies  of  the  micromecbani cs  of  failure  of  unidi rect i onal ly  reinforced 
composites  are  based  on  theoretical  analyses  of  the  stress  distribution  in 
damaged  composites  under  load.  Any  use  of  the  theory  of  elasticity  for  this 
purpose  turns  out  to  be  prohibitively  difficult,  and  therefore  approximate 
techniques  must  be  employed.  A  number  of  authors  have  used  shear  lag  theory 
or  some  modification  thereof  for  this  purpose  [1-5].  However  the  relation 
between  the  solutions  found  in  this  manner  and  experimental  results  has  been 
made  uncertain  by  the  fact  that  the  solutions  involve  a  parameter  the  values 
of  which  were  not  really  known.  This  parameter  is  the  effective  shear  stiff¬ 
ness  of  the  matrix,  which  serves  to  transfer  direct  stress  from  one  fiber  to 
another. 

Several  previous  papers  deal  with  various  aspects  of  this  question.  One 
gives  an  exact  solution  for  the  shear  interaction  between  two  infinitely  long 
rigid  fibers  immersed  in  an  infinite  elastic  medium  [6].  Another  gives  an 
approximate  solution  for  the  interaction  between  neighboring  pairs  of  fibers 
in  a  square  array  [7],  A  third  extends  this  solution  to  the  case  of  a  rec¬ 
tangular  array  [8]  . 

The  present  paper  considers  a  linear  array,  i.e.,  a  unidi rect ional ly 
reinforced  lamina  or  tape.  A  theoretical  solution  is  given  for  such  an  array 
that  is  valid  when  the  fiber  diameter  is  small  compared  to  fiber  separation 
and  fiber  separation  is  small  compared  to  matrix  thickness.  For  other 
geometries,  the  solution  is  obtained  experimentally  using  an  electric  analogue 
technique  [8],  These  findings  are  compared  to  assumptions  that  have  been  made 
concerning  the  shear  transfer  stiffness  by  previous  investigators. 
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Analytical  Solution  for  Linear  Array 

The  analysis  of  a  linear  array  can  be  carried  out  by  considering  the 

case  of  a  three-fiber  interaction  in  which  static  equilibrium  is  maintained  as 

in  Figure  1.  The  center  fiber  is  loaded  with  an  axial  line  force  2F  9  per 

o 

unit  length,  while  the  neighbors  are  each  loaded  with  F  ' ,  directed  in  the 
opposite  sense. 


The  shear  lag  equation  becomes 
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The  displacement  of  each  fiber  can  be  obtained  by  superimposing  the  effects  of 
all  of  the  applied  loads.  Employing  the  relationship  between  the  load  points 
and  the  field  point  in  [6], 
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Subtracting  (4)  from  (3),  we  obtain 

F  ' 

W  -  w  =-  [2  ln(d/r  )  +  ln(d/2r  ) 

a  b  2  nG  I  o  o 


Thus, 


2(wa  "  V  =  rfc  [3  ln(d/ro>  ~  ln2] 
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Comparing  (7)  with  (2), 


h  =  _ 2n _ 
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Experimental  Setup 

The  theory  underlying  the  electric  analogue  can  be  found  in  [8],  A 
linear  array  with  various  thicknesses  of  matrix  material  was  simulated  by 
placing  two  vertical  insulating  panels  with  height  exceeding  the  depth  of  the 
fluid  (electrolyte),  on  top  of  a  submerged  panel  with  holes  that  match  the 


size  of  the  conducting  rods  and  with  spacing  approximate  to  the  r  /d  under 


investigation*  The  insulating  panels  were  spaced  at  equal  distances  from  the 
center  lines  of  the  rods.  The  distance  was  chosen  to  represent  the  desired 


composite  thickness*  Thus,  the  fluid  between  the  vertical  panels  was  analo¬ 
gous  to  the  matrix  associated  with  that  particular  array. 


Linear  arrays  with  different  thicknesses  and  rod  spacing  were  investi¬ 
gated,  the  first  experiment  involving  three  rods  as  in  figure  2.  A  constant 
AC  potential  was  applied  across  rods  a  and  b  with  b  connected  to  c.  The 
current  through  the  power  supply  was  measured  in  each  case.  The  potential 
difference  between  a  and  b  was  also  measured.  The  governing  equation  can  be 
written  as 
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The  effect  of  the  presence  of  a  free  rod  adjacent  to  b  and  c  was  investigated 
by  a  second  experiment  as  in  Figure  3. 


Results 

The  results  of  the  experiments  for  evaluating  h/d  are  summarized  in  Fig¬ 
ures  4-6.  In  these  figures  h/d  is  plotted  against  r  / d  for  various  ratios  of 

fiber  radius  to  lamina  thickness,  r  / t.  We  do  not  have  a  theoretical  solution 

o 

except  for  r  /t— >0,  so  Figures  4-6  are  based  entirely  on  the  experimental 

data.  Extrapolation  beyond  the  data  points  is  aided  by  the  knowledge  that 

h/d  — >  0  for  r  /d  — »  0,  and  h/d  — >  •  when  r  /d  — >  0.5.  These  limiting  cases 
o  o 

are  discussed  in  Reference  7  and  8. 
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Figure  7  shows  the  results  for  h/d  for  a  linear  array  of  rods  in  an  open 

i 

j  tank  (i.e.,  the  vertical  panels  are  removed,  so  the  matrix  has  the  shape  of 

the  tank).  Also  shown  are  the  theoretical  results  for  an  infinite  tank.  We 

note  that,  as  expected,  the  theoretical  results  agree  very  closely  with  the 

|  experimental  data  obtained  in  a  large  but  finite  tank  for  small  values  of 

r  /d,  i.e.,  where  the  theoretical  solution  is  valid, 
o 

Another  limiting  case  is  obtained  by  noting  that  for  a  very  rectangular 
I  array  (a  >>  1)  the  nearest  fiber  interaction  should  be  almost  the  same  as 

though  the  adjacent  lamina  was  not  there.  Table  1  gives  the  h/d  obtained  for 
a  =  10.5  and  shows  that  the  results  are  substantially  the  same  as  those 
obtained  for  a  linear  array  in  an  open  tank. 

In  this  investigation,  most  of  the  data  were  obtained  using  three  rods 

as  in  Fig.  1.  A  linear  array  is  better  approximated  by  a  very  large  number  of 

rods,  but  a  very  large  number  of  rods  will  in  many  cases  exceed  the  capacity 

of  the  tank.  To  see  whether  additional  rods  are  really  needed,  the  effect  of 

adding  an  additional  rod  at  each  end  as  in  Fig.  3  was  investigated.  The  added 

rods  were  left  free  to  adopt  the  local  potential.  The  results  are  contained 

in  Table  2.  It  was  found  that  the  resistance  was  virtually  unaffected  by  the 

presence  of  the  additional  rods.  The  largest  change  was  a  couple  of  percent, 

and  it  occurred  when  r  / d  and  r  It  were  both  large.  If  either  was  small,  the 

o  o 

change  in  resistance  was  negligible.  We  conclude  that  the  results  obtained 
using  three  rods  are  sufficiently  accurate  for  most  engineering  purposes. 
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Some  investigators  who  have  employed  shear  lag  to  find  stress  distribu¬ 
tions  in  single  ply  composites  have  assumed  that  h/ d  is  equal  to  t/d, 
irrespective  of  the  value  of  r^/d.  Since  t/d  =  ( t/rQ) (ro/d) ,  a  plot  of  t/d 

vs,  r  /d  is  a  straight  line  of  slope  t/r  .  Such  lines  are  shown  in  dot-dash 
o  o 

format  for  r  /t  *  0.107,  0.144,  and  0,380.  We  see  that  the  assumption  that 
o 

h/d  =  t/d  is  not  bad  for  most  of  the  r  /d  range  when  r  /t  =  0.144,  but  is 

o  o 

rather  far  off  for  values  of  r^/t  differing  widely  from  this  value. 

This  paper  and  its  predecessors  [7,8]  have  attempted  to  find  the  shear 
interaction  to  be  used  under  the  simplifying  assumption  generally  employed  in 
shear  lag  theory  that  only  nearest-neighbor  fibers  interact.  The  accuracy  of 
this  basic  assumption  is  currently  unknown. 
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Table  2.  Resistance  due  to  3-fiber  interaction  and  5-fiber  interaction 
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FIG.  1  A  3-FIBER  INTERACTION  IN  WHICH  STATIC 
EQUILIBRIUM  IS  MAINTAINED 
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FIG.  2  SIMULATION  OF  A  3-FIBER  INTERACTION  BY 
AN  ELECTRIC  ANALOGUE 


FIG.  3  SIMULATION  OF  A  5-FIBER  INTERACTION  BY  AN  ELECTRIC 
ANALOGUE 
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